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We notice that VT symmetric non-Hermitian one dimensional simple Harmonic 
Oscillator under simultaneous transformation of co-ordinate and momentum with 
proper choice of positive oscillating frequency can reflect negative spectrum with well 
behaved wave function in real space .We also present a suitable comuter programme 
to realise the negative spectrum directly . 
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I. INTRODUCTION 

In Quantum Physics,understanding on two important things [1] :(i) well behaved 
wave function 

v k n (x —>■ oo) =>• 0 (1) 

and (ii) commutation relation 

[x,p]=i (2) 

is necessary .Here \l/ n is the nth state wave function of the Hamiltonian 


H I'Pn >= E n \V n > 


( 3 ) 


having energy eigenvalue E n .The best example to understand wave function is to 
visualise the problem of Harmonic Oscillator ,which is an interesting system in Quan¬ 
tum Physics .Its eigenvalue relation 



x = ix (7) 

prposed by Bender,Hook and Klevansky [2] . The Hamiltonian can be written as 

= + ( 8 ) 

whose energy eigenvalue is 

E n = —{n + -) (9) 

However the wave function represents 

** = (10) 

divergent nature in real space i.e 

T n (x —y oo) =>■ oo (11) 

Hence we believe that approach of BHK[2] is not a suitable representation in real 
space as one can not visualise the negative energy in real space in view of divergent 
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nature of wave function .In this context we notice that Fernandez[3] and Rath [4] 
independently discussed negative energy of Harmonic Oscillator under simultaneous 
transformation of co-ordinate x as [5] 


x —y 


and momentum p as [5] 


P 


x + iXp 

\/(l + /^) 

p + i/3x 


( 12 ) 


(13) 


\/(l + 

.It is worth to notice that both the authors have different opinion on visualising the 
negative spectrum .Rath [4] argues with positive frequency of oscillation wiorw 2 0 
,where as Fernandez [3] with negative frequency of oscillation u> <C 0 having divergent 
nature of wave function 


*l>o(x) =< a#o >= ~i/A e Ux2 ' 2 ( 14 ) 

There fore one is likely to believe that prposals of BHK[2] and Fernandez [3] even 
if suitable for complex space but can hardly be a plausible approach for real space . 
Hence in this paper we present a detail analysis of well behaved wave function in real 
space having negative energy of Harmonic Oscillator under simultaneous transforma¬ 
tion of momentum and co-ordinate .Further we also present a direct simple approach 
to realise the negative energy based on computational work . 

II. Invariance in commutation relation and transformed Hamiltonian in 
second quantisation relation. 


Now we consider the commutation relation 

(x + iXp) (p + i/3x) 


\J{ 1 + X/3) 


x,p\ = i 


(15) 


Now the new Hamiltonian with transformed x and p becomes non-Hermitian in nature 
and is 


(p + i/3x) 2 t (x + iXp ) 2 
2(1 + X/3) + 2(1 + X/3) 


2 



In order to solve the above Hamiltonian (Eqn. (7)), we use the second quantization 
formalism as 


(a + a + ) 

\/2u 


( 17 ) 


and 



(18) 


where the creation operator, a + and annhilation operator a satisfy the commuta¬ 
tion relation 


a, a + ] = 1 


(19) 


and u is an unknown parameter. The transformed Hamiltonian can be written as 


where 


and 


h = h d + h n 

(20) 

Id X2 V . . (l-/3 2 ) 1 (2a+a + l) 

'- [(1 AV+ a; 4(1 + A/3) 

(21) 

a 2 (a + ) 2 

(22) 

Hn ~ U 4(l + A f3) + ^4(1 + A (3) 

= [ u(l A 2 )+ (1_ ^ } 2(A + /3)\ 

OJ 

(23) 

[ cu(l A 2 )+ + 2(A + /3)] 

OJ 

(24) 


111(a) .Negative energy calculation (Case Study for U=0) 

Now we solve the the eigenvalue relation: 


H^ n (x) = e n ^ n (x) (25) 

using perturbation theory as follows. Here we express 

k 

= 4 0> + E 4 m) (26) 

m =1 
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The zeroth order energy e® satisfies the the following eigenvalue relation 


H d \ ip n >= H D \n >= e^\n > 


( 27 ) 


where i[)^p is the zeroth order wave function and e ^ is the mth order perturbation 


correction. 


and 


(0) = ( 2W + 1 ) [(1 _ A 2) + (1-£!)j 

4(1 + A/5) U ’ 0J J 


V e (m) = e (1) + e (2) + e (3) + 

/ v n n ' n ' n 1 

m= 1 


(28) 


(29) 


The energy correction terms will give zero contribution if the parameter is deter¬ 
mined from non-diagonal terms of Hn [5] 

Let the coefficient of a 2 is zero [5] i.e. 


U = [-w(l - A 2 ) + 


2\ , (1-/5 2 ) 


U 


+ 2(A + /5)] — 0 


which leads to 


In this case, 


oj = u i = 


(£- 1 ) 
(1 + A) 


$ ) = -(n + -) 

Now the perturbation correction term is 

H T r (« + ) 2 (A + fl 

w ^4(1 + A/5) 1 + A/5 ^ 

In this case one will notice that 

, , Jn(n - 1 ) 

<n|if„| n -2>= 


, + 12 


(30) 


(31) 


(32) 


(33) 


(34) 


< n — 2\H N \n >= 0 


(35) 
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Hence it is easy to note that all orders of energy corrections will be zero. Let us 
consider explicitly corrections up to third order using standard perturbation series 
given in literature [6-8], which can be written as 


4 1} =< ^n\H N \^ n >= 0 

J2) _ V' <^n\H N \ll) k ><tlj k \H N \tlJn> 

€ n ~ - , r ( 0 )_ (OK- 

_ < 1pn\H N \^ n+ 2 >< ^ n+2 \H N \^ n > 

(<t 0> - 4L) 

J3) _ <1pn\H N \'4’p><4’p\HN\4’q><i’q\HN\lpn> _ n 

— 2^p,q ~ U 

V c n ) v c n c q ) 


(36) 


(37) 


(3) _ < 'fpn | h/jy |'0n + 2 >< t^n+2 177jy |^ n +4 >< ^n+4|-fljv|^n > _ 


= 0 


(M (o) w (o) (o) \ w ( 38 ) 

ttn t n+2 ){tn t n+4y ) 

Here second order correction is zero due to < V’n|77jv|'0n+2 >= &n,n +4 and third 
order correction is zero due to < 'ip n+ 4 \H]y\'ip n >= 5 n+ ^ n+ 2 . Similarly one can notice 
all correction terms e ^ will be zero. Hence 

1, 


en = ei 0) = -(n+ ^ 


(39) 


\Yn >= ( 


yjii‘2 n n\ 


) 2 H n {y/u[x)e ^"2 


with 


< > = 1 

III.(b). Coresponding Wavefunction using Perturbation Theory 

Here we find the wavefunction as 

V J / n fc) = I i^n > +f\,l3^^-\4 ; n+2 > + (/a,^) 2 ^^^-|'0n+4 > + 

V( n + 6 ) ! | 


(h,pY 


48\/riT 


‘ lt^n+6 > +• 


where In its compact form one can write, 


= u-m^^]\l^^\^ + 2 k > W1 


k =0 


(1 + X/3) 


n\ 


(40) 


(41) 


(42) 


(43) 
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The normalization condition here can be written as[6-9] 


<?A|T«>=1 (44) 

So also the eigenvalue relation 

< il> n \H\V$ ] >= e n = -(n + i) (45) 

III.(c). Direct calculation of negative energy using MATLAB Here we 
would like to present a simple MATLAB programme [10,11] to demonstate directly 
with a view to show that above selection of positive frequency can lead directly 
negative energy . 


N = 100; s = 1; 

(46) 

n = 1 : N — 1; 

(47) 

m = ^[ n )\ 

(48) 

P = 2; 

(49) 

A = 0.5; 

(50) 

L ~ 1 + /3A’ 

(51) 

/5-1 

Z01= A + 1 ; 

(52) 

g 

-~i= -* ( diag(m , —1) + diag(m , 1); 

V(2»i) 

(53) 

v -j= — * (diag(vi, —1) — diag(m, 1); 

(54) 

R = p + l/fa; 

(55) 

y = x + iXp; 

(56) 

H = 0.5 * L(R 2 + F 2 ); 

(57) 

Air/ = sort(eig(H )); 

(58) 
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Eig{ 1 : 100 ); 


( 59 ) 


In Tablc-1 , we present first fifty -One eigenvalues of the transformed Harmonic 
Oscillator . 

Table -I : Negative Energy of transformed Harmonic Oscillator with 
positive frequency w± = 1,2, 3,4 


n 

A = 1, /3 = 3 

A = 2,/3 = 7 

A = 0.5,/3 = 5.5 

A = 0.8, /3 = 8.2 

0 

- 0.5 

- 0.5 

- 0.5 

- 0.5 

1 

-1.5 

-1.5 

-1.5 

- 1.5 

2 

- 2.5 

- 2.5 

-2.5 

- 2.5 

3 

- 3.5 

-3.5 

-3.5 

- 3.5 

4 

-4.5 

-4.5 

- 4.5 

- 4.5 

5 

-5.5 

-5.5 

-5.5 

- 5.5 

10 

- 10.5 

-10.5 

-10.5 

-10.5 

20 

-20.5 

-20.5 

-20.5 

- 20.5 

30 

-30.5 

- 30.5 

-30.5 

-30.5 

40 

-40.5 

-40.5 

-40.5 

- 40.5 

50 

- 50.5 

-50.5 

-50.5 

-50.5 


Here we notice eigenvalues are negative for different values of /3 and A 
with positive frequency of vibration . 

IV.(a). Negative energy calculation (Case Study for V=0) 

Let the coefficient of (a + )2 is zero [5] i.e. 

V = [-a;(l - A 2 ) + Q—P—L - 2(A + /3)] = 0 (60) 

id 


which leads to 


Ul = U>2 = 


(1 + P) 
(A — 1) 


Now the perturbation term becomes 


H n = U 


a 2 _ (A + /3) 

4(1 + A/3) 1 + A/3 


(61) 


(62) 
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In this case one will notice that 


< 4>n\HN\4>n+1 >~ U 


\J[{n + l)(n + 2)] 
4(1 + A 0) 


(63) 


< 4>n+2 1 H n | (j) n >— 0 


(64) 


Hence we have 


= < <l>n\H N \<j> n >= 0 


(65) 


4 2) = E 

k^n 


< (/>n\HN\4>k >< 0/c | Hpj | (j) n > < (j) n \H N \^) n -2 >< </>n-2|-^)v|0n > 


(ek L 


(0) ^(0)> 


(fn 


(0) (0) \ 
e n-2j 


= 0 

( 66 ) 


(3) | H N | (j)p > <4 (f)p | Hjv | (j)q > <4 (f)q | -Hjv 1^ _ 

£ ” =2j (4 0) -4" , )(4” , -4 0) ) = 


= o 


(67) 


P.<? 


or 


(3) ^ |-^iV | ^n—2 ^*4 0n—21^7V1 4^n— 4 0n—41-^JV1 p. ('CO') 

e?l = 716) (o) u (o) loTE = u ^ 

(^fcn — t n _ 2 )\tn — t n _ 4 J 

Here second order correction is zero due to < (j) n \H n\4 > u-2 >= &n,n- 2 and 
third order correction is zero due to < <p n - 4 \ 1 4> n >= <^n- 2 ,n- 4 - Similarly 
one can notice all correction terms ei m ) will be zero. Hence 


e„ = 4 0) = E^< = -{n + -) 


(69) 


and 


\K >= (-^L)iH n (^T 2 x) e -^ 


' y/n2 n n\ 

IV. (b). Coresponding Wavefunction using Perturbation Theory 


(70) 


Here we consider the wavefunction as 



*?’=!*• > >+ 

, 3 Vn! 


(/a,/?)" 1 - 

48^(n-6)! 

In its compact form, one can write 


|^n-6 > +• 


<f,(fe) _ 1 \fc/ ( A + 70 ^ I / 

/ y l f/ w-i | \ n\ ) r 77T . (pn—2k u )2 

^0 (! + A P) y(n — 2k)\2 k k\ 

Here we notice that for x —> 00 i.e. 

4> n (x OO) —> 0 

and 

—>■ 00 ) —* 0 

In this case, the normalization condition can be written as[6-9] 


< >= 1 


So also eigenvalue relation 


<K\H\^ >=E n = -(«+-) 


(71) 

(72) 

(73) 

(74) 

(75) 

(76) 


IV. (c). MATLAB programme to realise negative energy directly. Here 


one has to use above MATLAB programme directly with new positive 
frequency of vibration W 2 as 


P = 1; (77) 

A = 3; 

/3 + 1 

~ T > 
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(78) 

(79) 




In Table-II , we present first fifty -One eigenvalues of the transformed 
Harmonic Oscillator . 

Table -I : Negative Energy of transformed Harmonic Oscillator with 
positive frequency W\ = 1,2, 3,4 


n 

A = 3, /? = 1 

A = 2,/3 = 1 

A = 1.5, /3 — 0.5 

A = 1.5, /3 — 1 

0 

- 0.5 

- 0.5 

- 0.5 

- 0.5 

1 

-1.5 

-1.5 

-1.5 

- 1.5 

2 

- 2.5 

- 2.5 

-2.5 

- 2.5 

3 

- 3.5 

-3.5 

-3.5 

- 3.5 

4 

-4.5 

-4.5 

- 4.5 

- 4.5 

5 

-5.5 

-5.5 

-5.5 

- 5.5 

10 

- 10.5 

-10.5 

-10.5 

-10.5 

20 

-20.5 

-20.5 

-20.5 

- 20.5 

30 

-30.5 

- 30.5 

-30.5 

-30.5 

40 

-40.5 

-40.5 

-40.5 

- 40.5 

50 

- 50.5 

-50.5 

-50.5 

-50.5 


Here we notice eigenvalues are negative for different values of f3 and A 
with positive frequency of vibration . 

V. Conclusion 

In this paper, we suggest a simpler procedure for calculating energy 
levels and wave function of the non-Hermitian harmonic oscillator under 
simultaneous transfromation of co-ordinate and momentum using pertur¬ 
bation theory .Further direct study on negative energy using simple MAT- 
LAB programme support the analytical study using perturbation theory 
.more interesting part is the well behaved nature of wave function under 
positive frequency of vibration i.e 

3> 0, x —* oo) —> 0 (80) 
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<& n {w 2 S> 0 , x —* oo) — > 0 

which can hardly be seen in all previous work . 


( 81 ) 
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